A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 



ASTRONOMY 

AND 
ASTROPHYSICS 

February 1, 2008 



On star formation and chemical evolution 
in the Galactic disc 



L. Portinari and C. Chiosi 

^\ Department of Astronomy, University of Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy (portinari, chiosi@pd.astro.it) 

OS ■ 

0\ ; Received 10 June 1999 / Accepted 27 August 1999 



=3 : 

<: 
o ■ 

m : 



> 
m 

00 

o 

ON 

Oh 
I 

O 



Abstract. The abundance gradients and the radial gas 
profile of the Galactic disc are analysed by means of a 
model for the chemical evolution of galaxies. As one of 
the major uncertainties in models for galaxy evolution is 
the star formation (SF) process, various SF laws are con- 
sidered, to assess the response of model predictions to the 
different assumptions. Only some SF laws are successful in 
reproducing the metallicity gradient, and only if combined 
with a suitable infall timescale increasing outward (inside- 
out formation scenario). Still, it is difficult to reproduce 
at the same time also the observed gas distribution; we 
therefore suggest further improvements for the models. 

Key words: Galaxy: chemical evolution - Galaxy: abun- 
dance gradients - Galaxy: star formation 



1. Introduction 

Modelling and understanding the formation and evolu- 
tion of galaxies is a major issue of modern astrophysics, 
involving a variety of approaches: dynamical and hydro- 
dynamical simulations, chemical evolution models, stel- 
lar population synthesis techniques. Any galaxy model re- 
quires a recipe for the star formation (SF) process, which 
unfortunately is still poorly known, especially on the large 
scales relevant to galaxy evolution. A variety of SF laws 
has been suggested in literature especially for spiral galax- 
ies, where we can observe on-going SF on large scales. In 
chemical models for galactic discs, different prescriptions 
for the SF process are used by different authors. This is 
true also for models of the disc of the Milky Way, a system 
which we can observe and study in much closer detail than 
any other. It is often hard to compare models and conclu- 
sions by different authors, since it is not clear to what 
extent model results have a general validity or rather de- 
pend on differing assumptions for the SF law. 

The present work is addressed at testing different pre- 
scriptions for the SF law in the Galactic disc. To this aim 
we use a model for the chemical evolution of galaxies, able 



to handle different options for the SF law (Section ||). We 
review extant data on the abundance gradients and on the 
radial gas and star formation rate (SFR) profiles of the 
Galactic disc (Section ||). To these observational counter- 
parts we compare models with the various SF laws, trying 
to discriminate, if possible, a "best-fit" law (Section 
Summary and conclusions are finally drawn in Section^. 

A similar analysis was already performed by Prantzos 
& Aubert (1995), but an upgrade of their study is sug- 
gested since some basic observational results have mean- 
while changed, especially regarding the metallicity gradi- 
ent (§1). 

2. The chemical model and the SF laws 

The chemical evolution of the Galactic disc is simulated by 
means of the model of Portinari et al. (1998, hereinafter 
PCB98), to which we refer for a detailed description. Here 
we just overview the main features of the model. 

It is an open model with continuous infall, as suggested 
by the metallicity distribution of G-dwarfs in the Solar 
Vicinity (Lynden-Bell 1975) and by dynamical simulations 
of disc formation (e.g. Larson 1976). Open models also 
seem to be appropriate for galaxies more in general — 
for instance, see Bressan et al. (1994) for the analogous of 
the G-dwarf problem in elliptical galaxies. The infalling 
gas, assumed to be primordial, settles onto the disc at an 
exponentially decreasing rate with timescale T(r): 

o-j„/(r,i) = A{r)e~^ 

A{r) is fixed by the assigned final (present-day) surface 
density profile of the model disc, an exponential with 
scale-legth r^: 



Send offprint requests to: L. Portinari 



where ctq = (T{rQ,tG), total surface density at the Solar 
radius at the final Galactic age tc- 

In PCB98 the chemical model was applied only to the 
Solar Neighbourhood, while in the present paper several 
concentric cylindrical shells (rings) are calculated, span- 
ning the Galactic disc from 2 to 20 kpc. Each ring, typ- 
ically 2 kpc wide, consists of a homogenous mixture of 
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gas and stars which is assumed to evolve independently, 
neglecting any possible radial flows. The instantaneous re- 
cycling approximation (IRA) is relaxed in the calculation 
of stellar ejecta. The chemical evolution of the interstellar 
medium (ISM) in each ring r is described by the set of 
equations: 

^G,(r,<)= -X,{r,t)-^{r,t) + 

+ J^-■^{r,t~TM)R^{M)^{M)dM+ (1) 

where the Gi's are the surface gas densities of each chem- 
ical species i, normalized to the total surface density at 
the present age of the Galaxy a{r, to). The first term on 
the right-hand side represents the depletion of species i 
from the gaseous phase due to SF; the second term is the 
amount of species i ejected back to the ISM by dying stars; 
the third term is the contribution of the infalling gas. 

The various ingredients of Eq. (^ and its numerical 
solution are described in detail in PCB98. With respect 
to that model, the only difference is that here we adopt 
different options for the SFR, '^{r,t). We describe here 
below the various SF laws considered, together with their 
physical underpinnings. 

2.1. Schmidt law 

A few analytical prescriptions for the SF law in galactic 
discs are available in literature. The easiest physical law 
is Schmidt (1959) law: 

SFR cx 

where ag is the surface gas density. In spite of its crude- 
ness, Schmidt law still remains, after 40 years, the most 
popular SF law in galaxy models. 

Recent empirical determinations in spiral and star- 
burst galaxies indicate a "Schmidt exponent" k ~ 1.5 
(Kennicutt 1998). The observed correlation is quite good 
between the SFR and the total density of cold gas 
(III-I-H2), while no clear correlation apparently holds be- 
tween the SFR and the sole II2 component — at odds with 
expectations since stars form within molecular clouds. 
Therefore, we consider here a Schmidt law with exponent 
1.5, involving the overall gas density with no distinction 
between the atomic and the molecular phase. 

In the formalism of our model, we translate Schmidt 
law into a normalized SFR per unit surface density (the 
function in Eq. |l|) as: 



dG 



dt 



G-ir,t) (2) 



where we have introduced the constant normalization fac- 
tor cr(r0, <g)~*-''~"'^'' for the sake of expressing the SF efh- 
ciency 1/ in [t^^]. We adopt k = 1.5. 

Hereinafter, chemical models adopting the Schmidt 
law (|) will be labelled with S15. 



Other SF laws still assume that the SFR depends on 
the gas density, but also take into account additional phe- 
nomena affecting the SF process. As a result, the efficiency 
with which the available gas is turned to stars depends on 
the galactocentric distance r. Such laws can be written 
as "modified" Schmidt laws, with a proportionality coef- 
ficient, or efficiency, varying radially over the disc: 

SFR = iy{r) ag 

The detailed radial behaviour is related to the physical 
processes involved, as we discuss in the following cases. 



2.2. Spiral density waves trigger 

According to a theory dating back to Roberts (1969) and 
Shu et al. (1972), SF takes place in spiral arms due to 
gas compression in density waves. In this framework, Oort 
(1974) suggested that the SFR at a given galactocentric 
distance is determined by the rate at which the orbiting 
gas crosses spiral arms: 

SFR oc {n{r) - Qp) 

where ^{r) is the angular rotation velocity as a function 
of radius and flp is the angular velocity of the spiral pat- 
tern. In the scenario for disc galaxy formation within a 
pre-existing dark halo, and also according to viscous disc 
models (e.g. Saio & Yoshii 1990), the rotation curve is 
time-independent and keeps the same shape throughout 
the evolution. For a flat rotation curve typical of galac- 
tic discs ri(r) cx r~^, and since usually flp <C ^{r), we 
roughly get: 

SFR oc cr^ 

Such a SF law was adopted, for instance, by Wyse & Silk 
(1989) and by Prantzos & Silk (1998). 

A similar radial dependence of the SF efficiency was 
also obtained by Wang & Silk (1994), who treated the SF 
process as a consequence of gravitational instabilities in 
the disc. These are related to the local Toomre's stability 
parameter Q, itself related to the local epyciclic frequency 

Ke{f) OC i^{r) oc r^^ 

Normalizing the expression so that the efficiency co- 
efficient v is in units of [t~^], in our chemical model we 
express this "Oort SF law" as: 



*(r,t) 



a{r, to) 



G^{r,t) 



(3) 



_cr(ro,iG). 

Kennicutt (1998) noticed that his empirical relation 
between the current SFR and the gas surface density, 
yielding k ~ 1.5 if fitted with a Schmidt law, is equally 
well fitted with an Oort-type law with k = 1. Therefore, 
in the Oort SF law (|^) we will adopt an exponent n ~ 1.0. 
Such models will be labelled hereinafter with GIG. 

For the sake of comparison with other authors, we will 
though consider also the case of an Oort-type law with 
K = 1.5 (models 015, see § |!|). 
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2.3. Self-regulating gravitational settling + feed-back 



2.4- The Initial Mass Function 



Other prescriptions describe SF as a self-regulation pro- 
cess between the gravitational settling of the gas onto the 
plane of the disc and the energy feed-back from massive 
stars. The former phenomenon leads to compression and 
cooling of the gas layer, favouring its collapse into clouds 
and stars. Stellar winds and supernova explosions from 
massive stars provide a heating and turbulence source 
that supports the gas layer against further compression 
and cooling. The SFR results from the balance between 
these two opposing effects. Where the total surface mass 
density a is larger, gravitational settling and the related 
cooling are more efficient, and a larger SFR is required to 
reach the equilibrium feed-back. As a result, the SFR de- 
pends both on the gas surface density ag and on a. Since 
in spiral discs a displays an exponential profile, the SF 
efficiency inherits an exponential radial decline. Such a 
self-regulating SF process was first proposed by Talbot & 
Arnett (1975) as: 



SFR oc a{r)''-^a' 



(4) 



A similar formulation was suggested later by Dopita 
(1985). Self-regulating SFRs depending on the local 
surface mass density are also obtained in full chemo- 
dynamical models (Burkert et al. 1992). 

More recently, Dopita & Ryder (1994) discovered that 
an empirical relation holds in outer spiral discs between 
the Hq, emission, tracing on-going SF, and the surface 
I-band brightness, tracing the underlying old stellar com- 
ponent and thus the total surface density. This relation 
provides the observational counterpart to a SFR depend- 
ing on a. Dopita & Ryder derived a theoretical SF law 
(strongly reminiscent of the one by Talbot & Arnett): 

SFR oc a" cr™ (5) 

with n = 1/3, TO = 5/3, which fits very well with the ob- 
served relation. It can also account for the observed corre- 
lation between oxygen abundance and surface brightness 
in spiral discs (Ryder 1995). 

In the formalism of our chemical model, normalizing 
the SFR (^) so that the efficiency coefficient v is always 
in units of [t~^], we write: 



*(r,t) 



cr"(r,i) a"'-'^{r, to) 
cr{rQ,tG) 



n-^m— 1 



G™(r, t) 



(6) 



and we adopt n = 1/3, to = 5/3. Chemical models with 
this SF law will be labelled with DR. 

In PCB98 the SF law by Talbot & Arnett (Eq. |) 
was adopted. That SF law, with Schmidt-like exponent 
K = 1.5, naturally leads to similar predictions as models 
DR, since the two SF laws are analogous. Therefore, it's 
not worth in the following discussion to present also mod- 
els with Talbot & Arnett 's law separately. 



Besides the infall timescale r(r) and the SF efficiency ly of 
the various SF laws, another important model parameter 
is related to the Initial Mass Function (IMF). We adopt a 
power-law IMF: 

$(M) cx Af"^ 

where /i=1.35 for M < 2 M© as in Salpeter (1955) 
and fi=1.7 for M > 2 Mq following Scab (1986). Since 
the bulk of chemical enrichment is due to stars with 
M >1 Mq, a meaningful parameter is the fraction ( of 
the total stellar mass distributed in stars above 1 Mq. 



c 



100 



$(M) dM 



Fixing the "scaling fraction" C is equivalent to fixing the 
lower end M; of the IMF (see PCB98 and Chiosi & Mat- 
teucci 1982 for further details). In our models we always 
keep the IMF constant all over the disc and in time. 

3. Observational constraints 

In this section we review the available observational data 
on the abundance gradients and on the gas and SFR radial 
profiles in the Galactic disc. In the next section we will 
discuss how predictions from chemical models, adopting 
the different SF laws, compare to these observational data. 

3.1. Metallicity gradients 

Negative radial gradients of metallicity in spiral discs have 
long been established. Abundance data are usually de- 
rived from HII regions and from bright blue stars. These 
young objects mainly trace the present-day gradient for 
oxygen, and then for nitrogen, sulfur, neon and argon. For 
the Milky Way, the first extensive study of abundances in 
HII regions was presented by Shaver et al. (1983). Later 
studies based on optical spectra, in the inner and outer 
Galaxy, are by Fich & Silkey (1991) and Vilchez & Es- 
teban (1996). Abundances in HII regions have been also 
measured from FIR lines (Simpson et al. 1995, Afflerbach 
et al. 1997, Rudolph et al. 1997). Studies on HII regions 
generally agree on the existence of a Galactic oxygen gra- 
dient around -0.07 dex/kpc. Maciel & Koppen (1994) used 
Type II planetary nebulae as tracers of the present-day 
gradient, which they derive to be -0.06 dex/kpc, consis- 
tently with studies on HII regions. 

At odds with nebular studies, the stellar spectra of OB 
stars and associations seemed to suggest a very flat oxy- 
gen gradient (Gehren et al. 1985; Fitzsimmons et al. 1990, 
1992; Kaufer et al. 1994; Kilian-Montenbruck et al. 1994). 
The discrepancy disappears in the recent, extensive stud- 
ies by Smartt & RoUeston (1997) and by Gummersbach 
et al. (1998), both deriving a gradient of -0.07 dex/kpc. 
This revision of the oxygen gradient in our Galaxy, in fact. 
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Table 1. Observed abundance gradients of various elements. Ar is the galactocentric radial range covered by each 
resective study (in kpc). Gradients are expressed in dex/kpc. 



tracer reference 


Ar 


d[0/HJ dlJV/H| d{S/H\ 
dr dr dr 


HIT regions Shaver et al. (1983)'^' 
(optical) Fich & Silkey (1991) 

Vilchez & Esteban (1996) 


4-13 
12-18 
12-18 


-0.07±0.015 -0.09±0.015 -0.01±0.02 

— ~ — 
-0.036±0.02 -0.009±0.020 -0.041±0.020 


HIT regions Rudolph et al. (1997)'^' 
(FIR) Afflerbach et al. (1997) 


0-17 
0-12 


-0.079±0.009 -0.111±0.012 -0.079±-0.009 
-0.064±0.009 -0.072±0.006 -0.063±0.006 


OB stars Smartt & RoUeston (1997) 
Gummersbach et al. (1998) 


6-18 
5-14 


-0.07±0.01 — — 
-0.07±0.02 -0.08±0.02 — 


Type II PNcE Maciel & Koppen (1994) 


4-14 


-0.06±0.01 — -0.07±0.01 



Rcscalcd to tq — 8.5 kpc. 

Includes the data of Simpson ct al. (1995). 



suggested to us to revise as well the study by Prantzos & 
Aubert (1995) on the SF and chemical evolution of the 
Galactic disc. 

In spite of the overall agreement, it is still controver- 
sial whether the slope of the oxygen gradient is constant 
all over the disc or not. A flattening in the outer zones 
has been claimed by Vilchez & Esteban (1996) and Fich 
& Silkey (1991), while other studies do not suggest such 
a trend (Smartt & Rolleston 1997; Rudolph et al. 1997; 
AfHerbach et al. 1997). In the overall, considering all the 
different data and tracers together (Fig. we consider 
that a single-slope gradient is a good enough description, 
and take 

= _o.06 - -0.07dexkpc-i 

dr 

as our estimate of the oxygen gradient throvighout the 
Galactic disc. 

A similar controversy also exists for the [N/H] gradient 
(Fig H top panel) which flattens out according to optical 
studies (Shaver ct al. 1983; Fich & Silkey 1991; Vilchez & 
Esteban 1996), while it remains roughly constant accord- 
ing to FIR measurements (Simpson et al. 1995; Rudolph 
et al. 1997; Afflerbach et al. 1997). The problem seems 
to arise from the very different, by 0.3-0.4 dex, nitrogen 
abundance determinations around r ~ Tq in the two sets 
of analysis (Rudolph et al. 1997). 

The very flat gradient of sulfur found by Shaver 
et al. (1983) is not confirmed in later studies; the sul- 
fur gradient (Fig. ^, middle panel) seems rather to fol- 
low the oxygen gradient, or equivalently the [S/0] ra- 
tio seems to be constant with radius, as expected for an 
abundance ratio of two primary elements produced by the 
same source (type II SNas). There seems to be no longer 
need to account for a gradient in the [S/0] ratio by invok- 
ing slightly different progenitors for the two elements, like 
SNae from massive stars of different mass ranges (Mat- 
teucci & Frangois 1989). 

The abundance gradients as determined in the above 
mentioned studies are listed in Table ^ the corresponding 
data are displayed in Fig. |l|. For the sake of completeness. 



from the above mentioned studies we plot also the data 
concerning the Galactic Centre. However, we will not con- 
sider these data as a constraint for chemical models of the 
disc, since the stellar population of the Galactic Centre is 
likely to have evolved on its own, rather than represent an 
extension of the disc population to r ^ 0. Our disc models 
will consider only regions down to about r = 2 kpc, as fur- 
ther inward the Bulge becomes the dominating Galactic 
component. 

The radial gradient of [Fe/H] is mainly traced by open 
clusters. An up-to-date determination based on a homo- 
geneous sample of open clusters can be found in Carraro 
et al. (1998), who find 

^[^M^_0.07--0.09dexkpc-i 
dr 

in good agreement with previous results (Janes 1979; Friel 
1995). 

3.2. Gas and SFR radial profiles 

The radial profile of the cold gas distribution (sum of 
HI and H2) is taken from the review by Dame (1993; 
Fig. ^, top panel). For r > 7 kpc, the profile is shallow 
with a plateau or a mild decline outward, while around 
4-7 kpc a pronounced peak is due to the molecular ring. 
For r < 3 kpc, there is a rapid drop in the gas density. 

Estimates for the SFR over the disc come from the dis- 
tribution of SN remnants (Guibert et al. 1978), of pulsars 
(Lyne et al. 1985), of Lyman-continuum photons (Giisten 
& Mezger 1982) and of molecular gas (Rana 1991 and 
references therein). Following Lacey & Fall (1985), since 
these observables cannot directly yield the absolute SFR 
without assumptions on the IMF etc., we just normalize 
them to their value at the Solar ring, and consider them 
as tracers of the radial profile SFR(r)/SFR(r0) (Fig. |, 
bottom panel). Unsurprisingly, the observed SFR profile 
follows that of the gas distribution, with a peak around 
r = 4 kpc and a milder decrease outward of the Solar ring. 
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Table 2. Parameter values and resulting metallicity gradients of the models. 



QT? lour 


u 


/■ 
t, 


Td [kpc] 


model 




d^ [dex kpc ] 


S15 


0.35 


0.2 


4 


S15a 
S15b 


3 (uniform) 

3e '■d 

(from r ~ 0.7 at r = 2.5 to r ~ 50 at r = 20) 


-0.03 
-0.055 


OlO 


0.19 


0.2 


4 


OlOa 
OlOb 


3 (uniform) 

T-TQ 

3 e ""d 

(^ironi T u. ( at r — z.o to t ou at r — zuj 


-0.03 
-0.06 


DR 


0.45 


0.2 


4 


DRa 
DRb 

DRc 


3 (uniform) 

3 (r > ro) 
~ 0.1 + '■"^■^ (r < rp)) 
(from r ~ 0.1 at r = 2.5 to r = 3 at r > Tq) 

3 (r > 4.5 kpc) 
~ 1 + (r- - 2.5) (r < 4.5 kpc) 
(from r ~ 1 at r = 2.5 to r = 3 at r > 4.5) 


-0.07 (r > re) 
flat (r < r©) 

-0.07 
-0.06 


015 


0.35 


0.2 


4 


Uioa 
015b 


o (umtorm j 
~ 3^ 

(from T ~ 1 at r = 2.5 to r ~ 7 at r = 20) 


-O.Uo > tq) 
flat (r < r©) 

-0.065 


DRr9 


0.95 


0.25 


4 


DRrQa 
DRr9b 


9 (uniform) 

9 e 3 fcp<: (r < 10 kpc) 
~ 50 (r > 10 kpc) 


-0.06 (r > 12 kpc) 
flat (r < 12 kpc) 

-0.06 


S15rd35 


0.35 


0.2 


3.5 


S15rd35a 
S15rd35b 


3 (uniform) 

3e '■d 

(from T ~ 0.7 at r = 2.5 to r ~ 50 at r = 20) 


-0.035 
-0.06 


O10rd35 


0.19 


0.2 


3.5 


O10rd35a 
O10rd35b 


3 (uniform) 

3e --d 

(from r ~ 0.7 at r = 2.5 to r ~ 50 at r = 20) 


-0.03 
-0.06 


DRrd35 


0.45 


0.2 


3.5 


DRrd35a 
Drd35Rb 


3 (uniform) 

3 (r > 6.5) 
~ 1 + (r < 6.5) 
(from r ~ 0.1 at r = 2.5 to r = 3 at r > tq) 


-0.075 (r > re) 
flat (r < r©) 

-0.075 



4. Model results 

In this section we compare model results for the various 
SF laws to the observed radial profile of the disc. 

The models are calibrated as follows. The age of the 
Galaxy is tG=15 Gyr. The final density profile of the disc 
air^tc) is exponential with scale-length and with a 



surface mass density at the Solar ring (r© = 8.5 kpc) of 
cr(r0, ic) = 50 Mq pc^^ (PCB98 and references therein). 
In most models, = 4 kpc and the infall timescale at the 
Solar ring is tItq) = 3 Gyr; different values for the infall 
timescale at the Solar ring and for the scale-length are 
discussed in § 4.7 and § 4.8. The adopted IMF is from 
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Fig. 1. Predicted radial oxygen gradient from the various models (see text) compared to observational data. The 
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Scalo (1986), uniform over the disc and constant in time, 
with an "amphtude" factor for SNae la, A = 0.07 (see § ^.41 
and PCB98 for further details). For each distinct SF law, 
the SF efficiency v and the "IMF scaling fraction" ( are 
calibrated so that at the Solar ring the model recovers 
the observed present-day surface gas density and oxygen 
abundance in the ISM (8-9 Mq pc"^ and [0/H]~ 8.7 dex, 
respectively). With respect to the latter constraint, no- 
tice that the local ISM is oxygen-poor with respect to the 
Sun (e.g. Meyer et al. 1998, and Fig. ||), at odds with ex- 
pectations from usual chemical models where metallicity 
always increases and should be higher in the present-day 
ISM than in the Sun, formed about 5 Gyr ago. Different 
explanations have been put forward to account for this 
effect (dust depletion, orbital diffusion, SN II pollution 
at Sun formation, recent metal-poor infall...), but all of 
them are rather ad hoc or uneasy to model properly. Here 
we will not deal with this problem; following Prantzos & 
Aubert (1995) we just concentrate on the radial metal- 
licity gradient, irrespectively of what the true zero-point 
abundance at the Solar ring actually is. Since the oxygen 
gradient is measured in young objects, we prefer to com- 
pare alike with alike and calibrate the model onto the ^^O 
abundance of the local ISM, rather than onto the solar 
one. Needless to say that models so calibrated will not 
reproduce the local age-metallicity relation and the solar 
metallicity; but the particular choice of the zero-point will 
not affect the radial behaviour of the models. 

Once we have so fixed the same Solar point for all the 
models, we analyse the radial behaviour in the various 
cases. The calibrated values for v and C, for the various 
models presented are given in Table |2| together with model 
predictions for the oxygen gradient. 

First of all we consider the constraint provided by the 
oxygen gradient as the best tracer of the overall metallicity 
(Wheeler et al. 1989), since it is the most abundant metal 
and the best sampled element in tracers of the present- 
day gradient (§||). Besides, its production and theoretical 
yields are the best understood of all the metals (Prant- 
zos 1998), so the modelling of the abundance evolution of 
oxygen is quite reliable. 



4-.1. Models with Schmidt law (S15) 

When a Schmidt law is adopted (Eq. ^ with k = 1.5), if 
no radial variation of model parameters is assumed, the 
predicted gradient is far too flat (model S15a, Fig. |l|). This 
is a well known result in literature (e.g. Lacey & Fall 1985; 
Edmunds & Greenhow 1995). Some further assumption is 
needed in models S15 to get a steeper oxygen gradient. 

A typical prescription is to assume an "inside-out" for- 
mation scenario for the Galactic disc, namely that the in- 
fall timescale r(r) is an increasing function of r (Chiosi & 
Matteucci 1980; Lacey & Fall 1983; Matteucci & Frangois 
1989). This can formally allow the model to predict a 



gradient close to the observed one (model SlSb, Fig. p, 
but the required variation of the infall timescale from the 
centre to the outskirts of the Galaxy is too large to be 
physically acceptable: in model S15b we need timescales 
much longer than the age of the Universe (r ~ 50 Gyr at 
r = 20 kpc, see Table H). This would imply that the accre- 
tion process is still largely on-going; most of the Galaxy 
mass on the outskirts would be still accreting onto the disc 
nowadays. Such high accretion rates in spiral galaxies at 
the present day are not observed in the X-rays (Berlin & 
Lin 1996 and references therein), meaning that the accre- 
tion process should be mostly over by now. 

We conclude that a simple Schmidt law is hardly suited 
at reproducing the observed metallicity gradients, at least 
within this kind of chemical models. 

4.2. Models with Oort-type SF law (OlO) 

In the case of the SF law by Oort (Eq. |^ with n ~ 1), 
model OlOa with a uniform infall timescale t — 3 Gyr 
still predicts too flat a metallicity gradient. To recover the 
observed gradient in the inside-out formation scenario, as 
in the case of models 515 the required radial variation 
of T becomes too extreme to be acceptable (model OlOa, 
Table H). This kind of SF law also proves unable to re- 
produce the observed oxygen gradient, at least within the 
framework of these chemical models. 

4.3. Models with Dopita & Ryder's SF law (DR) 

Model DRa is run with a uniform infall timescale r=3 Gyr 
all over the disc. The predicted oxygen gradient is in good 
agreement with the data in the outer regions of the disc, 
while it flattens considerably in the inner regions (Fig. |l|). 




5 10 15 20 

R(kpc) 



Fig. 2. Evolution of the oxygen gradient in model DRa. 
The gradient is plot at age steps of 1 Gyr, from 1 to 15 
Gyrs. 
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It is worth commenting a bit further on this effect. A 
hnear gradient in [0/H] = 12+log(0/H) impUes an expo- 
nential decline of the metallicity Z. Therefore, a SF effi- 
ciency declining exponentially in radius, such as that by 
Dopita & Ryder, looks appeaUng to account for the metal- 
licity gradient. In fact, simple closed models with such a 
SF law easily predict an exponential radial decrease of Z 
(Edmunds & Greenhow 1995). But in more realistic mod- 
els this prediction falls, and the metallicity in the inner 
regions tends to settle onto some equilibrium value rather 
than follow the espected exponential profile. As noticed by 
various authors (Gotz & Koppen 1992; Ferrini et al. 1994; 
Prantzos & Aubert 1995; Ryder 1995; MoUa et al. 1996, 
1997), when infall is included and the IRA is relaxed, the 
model reaches a stage in which the metal enrichment pro- 
vided by the ejecta of massive stars is significantly diluted 
by the primordial infalling gas and by the gas expelled 
by long-lived, low mass stars which formed at early times 
from relatively metal-poor gas. At this stage, metallicity 
abandons the one-to-one increase with SF processing and 
tends to settle on some equilibrium value. This effect is 
best seen in the inner regions of the Galaxy, where the SF 
process is more efficient: at the present age, most of the 
gas has been eaten up already and SF has so much re- 
duced as to be compensated by the gas shed by the many 
low-mass stars which had formed in the early, very active 
phases. As a consequence, in realistic models the radial 
metallicity profile tends to flatten in the inner regions. 
Fig. ^ shows how this effect develops in time: the gradi- 
ent is linear (in [O/H]) at early ages, while for t > 5 Gyr 
the metal enrichment slows down in the central parts due 
to the dilution effects of infall and of delayed return from 
low-mass stars. 

Therefore, even in the case of an exponentially varying 
SF efficiency as in models DR, some further assumption 
is needed to steepen the abundance profile in the inner 
regions. In the inside-out formation scenario, we take the 
infall process to be faster in the inner regions, while in 
the outer regions (r > tq) we can just keep a uniform 
r = 3 Gyr (model DRb, Fig. |l|, dashed line). 

This kind of SF law, combined with an inside-out for- 
mation scenario, can reproduce the observed oxygen gra- 
dient within acceptable infall timescales over the whole 
disc. 



4.4- intermediate case: models 015 

Though with an Oort-type SF law the preferred empiri- 
cal exponent is suggested to be k = 1 (Kennicutt 1998), 
some authors adopt k = 1.5, which may be still compati- 
ble with Kennicutt 's observations provided the whole disc 
model is considered altogether (Boissier & Prantzos 1999; 
Prantzos & Silk 1998). For the sake of completeness, we 
briefly discuss this SF law as well (Eq. |3| with k = 1.5). 

The behaviour of models 015 is qualitatively interme- 
diate between models S15 or OlO and models DR. Model 



015a, adopting a uniform infall timescale r = 3 Gyr all 
over the disc, predicts a metallicity gradient which is 
steeper than in models S15a or OlOa, but still shallower 
than the observed one, and shallower than in model DRa. 
The same flattening trend in the inner regions as in model 
DRa is also observed. With this SF law, model 015b is 
tuned to reproduce the observed gradient by means of a ra- 
dial increase of the infall timescale r, in the inside-out for- 
mation scenario (cfr. Table At odds with model DRb, 
here r(r) needs to increase over the whole radial range, 
otherwise the predicted gradient is not steep enough in 
the outer regions; but, unlike models S15b or OlOb, one 
does not need to resort to unplausible timescales in the 
outskirts to recover the observed gradient. 

Therefore, models 015 are acceptable toy-models (e.g. 
Prantzos & Silk 1998). We will not discuss them any fur- 
ther though, since we feel that an Oort-type SF law should 
be combined with the empirical exponent k = 1 rather 
than 1.5 (Kennicutt 1998). In all what follows, the reader 
might simply keep in mind that models 015 would qualita- 
tively display an intermediate behaviour between models 
S15/O10 and models DR. 
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Fig. 3. Predicted nitrogen, sulfur and iron gradients from 
models indicated in the legend on bottom left (see text) 
compared to observational data. 
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4-.5. Modelling the gradients of other elements. 

We have so far discussed the gradient of oxygen as the 
best tracer of the overall metallicity. With respect to the 
radial gradients of the other elements (§ ||), all the mod- 
els presented above which reproduce the observed oxygen 
gradient (S15b, OlOb, DRb) are also successful in predict- 
ing the other gradients (Fig. ||). This simply means that 
the adopted nucleosynthetic yields (from PCB98) are re- 
liable. In fact, once a model is calibrated in i/, ( and t so 
as to give the right oxygen abundance and gradient, the 
predicted behaviour of other elements merely depends on 
the relative abundance ratios implied by the yields. 

Fig. 1^ just shows some problem with the zero-point of 
the sulfur gradient (middle panel): the slope of the gra- 
dient is correct but the overall abundance is too high. 
This can be ascribed to the uncertainties in the yields of 
^^S from SNae II. The mismatch between the predicted 
and observed sulfur abundance goes here in the opposite 
direction with respect to evidence from stars in the So- 
lar Neighbourhood, which rather indicates that theoreti- 
cal yields for '^^S are too low (PCB98). But the data on 
the ^^S abundance in low- metallicity local stars are quite 
few and old; further and improved observational data are 
needed to shed light on the problem of the yields of ^^S. 

On the base of the constraint provided by the metallicity 
gradient, we can consider DR as viable models, while we 
tend to exclude models S15 and OlO, and the correspond- 
ing SF laws, since they would require unacceptably long 
accretion timescales in the outskirts of the Galaxy. The 
situation may though change if radial inflows of gas are 
allowed for (e.g. Lacey & Fall 1985; Gotz & Koppen 1992; 
Portinari & Chiosi 1999). 

4-. 6. Modelling the gas and SF profiles 

All of the models which are tuned, with more or less plau- 
sible model parameters, to give successful predictions for 
the metallicity gradient (S15b, OlOb, DRb) fall short at re- 
producing the observed radial profiles of gas density and 
SFR (Fig. ^). In particular, they cannot account for the 
peak in the gas density corresponding to the molecular 
ring around 4 kpc, and the corresponding peak in SF ac- 
tivity. Model S15 also seems to predict too high a SFR 
outward of the Solar ring, in agreement with the remark 
by Prantzos & Silk (1998) that the present-day SFR along 
the disc does not resemble any simple power law of the gas 
density. This latter feature favours, in general, SF laws 
with a radially varying efficiency (as in models DRb and 
OlOb) with respect to simple Schmidt laws. 

In models S15b, OlOb, DRb the radial dependence of 
the infall timescale was tuned so as to optimize the slope 
of the predicted metalhcity gradient (Fig. |^); these models 
tend to yield a maximum in the gas distribution around 
r0 (Fig. ^). One can also adjust the infall timescale so 



Dame (1993) DRb 
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Fig. 4. Observed radial gas and (normalized) SFR distri- 
bution, together with predictions from the models in the 
legend on top right. 

as to shift the maximum in gas density and SFR around 
4-5 kpc, while keeping somewhat compatible with the ob- 
served metallicity gradient (e.g. model DRc in Fig.|5|); still, 
the sharp peak of the molecular ring is never recovered (see 
also Prantzos & Silk 1998). 

It seems to be impossible for these models to reproduce 
both the metallicity gradient and the radial gas distribution 
at the same time. This is easily interpreted in the frame- 
work of "static" chemical models, or models consisting of 
isolated rings, like ours. To produce a higher metallicity 
in the inner regions, SF must have been more effective 
there, and consequently the available gas must have been 
more efficiently processed, consumed and locked in long- 
lived stars or remnants. In the models, the present-day 
gas density decreases where the metallicity is higher, at 
odds with observations. 

A possible way out, in the framework of simple static 
models, might be to assume that the IMF is skewed to- 
ward high masses in the inner regions (higher C. or flat- 
ter power-law exponent). A higher percentage of massive 
stars provides in fact a more efficient chemical enrichment, 
while less gas remains locked in low-mass stars and the 
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Dame (1993) 
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Fig. 5. Model DRc, tuned so as to predict a maximum in 
the gas distribution around 4-5 kpc, still cannot reproduce 
the observed peak (see text) 



gas density can keep high. But the issue of systematic 
variations in the IMF is still very controversial, both on 
empirical and on theoretical grounds (Scalo 1998a, b). 

A less "exotic" possiblity than a strongly varying IMF 
along the disc is that the formation of the molecular ring 
involves some kinematical or dynamical process, not ac- 
counted for in simple static chemical models. As is the 
case for rings in outer barred spirals, it is likely that the 
molecular ring in our Galaxy has formed because of gas ac- 
cumulation in correspondancc to th(^ Lindblad resonances 
of a barred potential. Models for the Galactic bar set its 
CO rotation between 2.5 and 3.5 kpc, and its Outer Lind- 
blad Resonance between 4.5 and 6 kpc (e.g. Gerhard 1996 
and references therein), which suggests a correlation be- 
twcion the molecular ring and the dynamical effects of the 
bar. In order to reproduce these effects, chemical models 
should allow for Bar-induced radial drifts of gas. Hereon 
the need to develop chemical models with radial flows of 
gas, which will be the subject of a second paper (Portinari 
& Chiosi 1999). 



O Gummersbach el al. (1998) 
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Fig. 6. Predicted oxygen gradient and gas distribution 
from model DRrQa adopting a uniform infall timcscalc 
T = 9 Gyr (solid line), compared to model DRa adopt- 
ing r = 3 Gyr (dotted line) and to the observational data. 
Model DRrQb, adopting r = 9 Gyr at the Solar ring, is 
calibrated with r increasing with radius so as to reproduce 
the observed gradient (see text). 

J^.l. The effects of a longer infall timescale 

All the models presented so far adopt an infall timescale 
of 3 Gyr at the Solar ring, as the fixed point for the ra- 
dial variation of r in the inside-out formation scenario. 
The accretion timcscalc for open chemical models of the 
Solar Neighbourhood is usually fixed by the metallicity 
distribution of local G-dwarfs. For many years, the refer- 
ence datasct was that by Pagel & Patchett (1975), later 
revised by Pagel (1989) and by Sommer-Larsen (1991a). 
This distribution was nicely fitted with infall timescales 
of 3-4 Gyr (e.g. Chiosi 1980; Matteucci & Prangois 1989; 
Sommer-Larsen 1991a). 

The latest datasets by Wyse & Gilmore (1995) and 
Rocha-Pinto & Maciel (1996) display a narrower distri- 
bution of the stars in metallicity, with a prominent peak 
around [Fe/H]~ —0.2, and therefore favour much longer 
accretion timescales, up to 8-9 Gyr (Chiappini et al. 1997; 
PCB98). Such timescales may be uncomfortably long from 
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Fig. 7. Top panel: surface mass density profile (solid line) 
and corresponding stellar density profile (dotted lines) for 
models with = 4 kpc. Bottom panel: the same plot for 
models with = 3.5 kpc. 



the dynamical point of view: in dynamical simulations the 
gas tends to settle onto the equatorial plane within the 
first 3-6 Gyrs (e.g. Burkert et al. 1992; Sommer-Larsen 
1991b; although Samland et al. 1997 suggest longer disc 
formation timescales). 

Even more problematic here, if the adopted infall 
timescale is so long, it is hard for the chemical model to re- 
produce the observed metallicity gradient. Let's consider, 
among the various SF laws studied here, the one which in- 
trinsically gives the steepest metallicity gradient, namely 
Dopita & Ryder's law. Model DRr9a adopts a uniform 
infall timescale r = 9 Gyr, rather than 3 Gyr (Table 
The resulting metallicity gradient is quite shallow (Fig.^, 
solid line), much shallower than the gradient predicted by 
the analogous model DRa with t — i Gyr (dotted line). In 
fact, when the infall timescale is very long, the dilution ef- 



fect (§ 4.3) is stronger and an equilibrium value for metal- 



licity is reached over a larger range of radii, leading to an 
overall flat gradient. To reproduce the observed gradient 
within the inside-out scenario (model DRr9b, dashed line) 
the assumed infall timescale for r > 10 kpc needs to be of 



Fig. 8. Predicted oxygen gradient and gas and SFR pro- 
files from model S15rd35a adopting a disc scale length 
Td — 3.5 kpc (solid line), compared to model S15a adopt- 
ing rd — 3.5 kpc (dotted line) and to the observational 
data. Model S15rrf35b is calibrated with r increasing with 
radius so as to reproduce the observed gradient (see text). 



^50 Gyr (Table [^), implying too large infall rates at the 
present day. With other SF laws, things get even worse. 

Radial inflows of gas in the disc can improve on this 
point, since they both affect the predicted G-dwarf dis- 
tribution contributing to solve the G-dwarf problem, and 
favour the formation of the metallicity gradient (Lacey & 
FaU 1985; Clarke 1991; Gotz & Koppen 1992). This wiU be 
discussed in more detail in a forthcoming paper (Portinari 
et al. 1999, in preparation). 

Finally, Fig. ^ shows how it is still impossible to re- 
produce at the same time the metallicity gradient and the 
gas peak around 4 kpc, also with models DRrQa.b. Namely, 
the conclusions drawn in 



4.f: remain valid irrespectively 



of the adopted T{rQ). 



4-.8. The effects of a shorter disc scale-length 

All the models presented so far adopt an exponential scale 
length for mass density = 4 kpc, quite typical of chemi- 
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Fig. 9. Same as Fig. pL but for models OlO. 
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Fig. 10. Same as Fig. ^ and ^, but for models DR. 



cal evolution models (e.g. Sommer-Larsen 1991b, Koppen 
1994, Rana 1991) and consistent with early determination 
of the optical scale length of the Galactic disc (e.g. Bin- 
ney & Tremaine 1987 and references therein; Kuijken & 
Gilmore 1989). Recent observations, however, tend to sug- 
gest a shorter scale length for disc stars, < 3 kpc (Sackett 
1997 and references therein, Vallenari et al. 1999 and ref- 
erences therein), comparable to that of other Sbc galax- 
ies. With rd ~ 4: kpc, our chemical models yield a final 
(present-day) stellar density profile with scale length be- 
tween 3 and 3.5 kpc (Fig. top panel), still compatible 
but somewhat long with respect to the latest observations. 
Therefore, in this section we explore models with a shorter 
mass density scale-length, = 3.5 kpc, corresponding to 
a final stellar density scale length around 2.5-3 kpc (Fig. 0, 
bottom panel). 

Models S15rrf35a, OlOr^BSa and DRr^SSa correspond 
exactly to models S15a, OlOa and DRa respectively, a part 
from the adopted scale length (see Table ^). The shorter 
scale length tends to favour the formation of the metallic- 
ity gradient and to predict a bit steeper decline of the gas 
density in the outer parts (compare solid to dotted lines in 
Fig. H to |lo|) . The effect is though minor on models S 15 and 
OlO, so that extreme infall timescales still need to be in- 



voked to reproduce the metallicity gradient in the inside- 
out scenario (models S15rc;35b and O10rc;35b in Figs. || 
and ^ see also Table The effect of the steeper disc pro- 
file is more evident in models DR, where the SF efficiency 
scales directly with surface density (§ 2.3); nontheless, to 
reproduce the observed metallicity gradient shorter infall 
timescales need to be assumed in the inner regions, in 
accordance with the inside-out scenario. The problem of 
reproducing the molecular ring is still present as well. 

In the overall, adopting a shorter scale-length (within 
observational limits) does not influence the bulk of the 
conclusions we drew in the previous sections. 

5. Summary and conclusion 

From our analysis of chemical models for the Galactic disc 
with different SF laws, it is evident that none of the SF 
laws considered is able, by itself, to reproduce the observed 
metallicity gradient throughout the whole extent of the 
disc. Some further "dynamical" assumption is needed in 
any case, such as an inside-out formation scenario (Lar- 
son 1976; Chiosi & Matteucci 1980; Lacey & Fall 1983; 
Matteucci & Francois 1989). 

Even SF laws with an exponentially decreasing effi- 
ciency (e.g. Dopita & Ryder's law), which in simple closed 
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models would predict a logarithmic metallicity gradient 
like the observed one, display a more complex behaviour 
in detailed open models with no IRA. In particular, due 
to the late dilution of the enriched ejecta of massive stars 
with infalling gas and returned gas from long-lived stars, 
the metallicity gradient tends to flatten in the inner re- 
gions with respect to the outer slope. To recover the ob- 
served slope, infall timescales need to be shorter in the in- 
ner regions, to reduce the effects of dilution at late times. 
This shows how relaxing the IRA has non-negligible effects 
on chemical models, even in the case of oxygen which is 
produced on very short timescales and for which it is of- 
ten assumed that neglection of finite stellar lifetimes is 
a viable approximation. Though this might hold for the 
massive stars producing oxygen itself, the gas shed on de- 
layed return from low-mass stars ultimately affects also 
the evolution of the oxygen abundance. 

A Schmidt SF law with k = 1.5 or an Oort-type SF 
law with K = 1.0, though nicely fitting some empirical 
relations (Kennicutt 1998), are unable to account for the 
observed metallicity gradients even in the inside-out pic- 
ture, because the required radial variation of the infall 
timescale T(r) implies accretion timescales in the outer 
regions much longer than a Hubble time. These SF laws 
are therefore excluded, within the framework of simple 
"static" chemical models. These conclusions may change 
in the presence of radial flows of gas (Portinari & Chiosi 
1999). 

We verified that these results are valid in general, ir- 
respectively on the adopted scale length of disc surface 
density (at least within observational estimates) and on 
the adopted infall timescale at the Solar ring. If the latter 
is very long, as maybe suggested from recent results on the 
"G-dwarf problem" , even SF laws of the Dopita & Ryder 
type hardly reproduce the metallicity gradient within the 
inside-out scenario. 

Finally, with these "static" chemical models it turns 
out to be very difficult to reproduce at the same time 
both the metallicity gradient and the gas distribution, in 
particular the molecular ring (unless one resorts to exotic 
possibilities like a varying IMF along the disc). Since the 
molecular ring is likely to be due to the dynamical in- 
fluence of the Galactic bar, chemical models allowing for 
radial drifts and accumulation of gas might be needed to 
mimick these effects and reproduce the corrects gas dis- 
tribution. Our next step will be to improve the chemical 
model including radial flows of gas (Portinari & Chiosi 
1999). 
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